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In the situation of two electro-statically coupled conductors a current in one conductor may 
induce a current in the other one. We will study this phenomenon, called Coulomb drag, in the 
Hubbard chain where the two "conductors" are given by fermions with different spin orientation. 
With the aid of a Monte Carlo (MC) approach which we presented in a recent paper we calculate the 
Transconductance in different variants of the Hubbard chain (with/without impurity and additional 
[long-ranged] interactions) for different fillings. 

PACS numbers: 75.30.Gw,75.10. Jm,78.30.-j 
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I. INTRODUCTION 



The Coulomb-drag effect describes how two conduc- 
tors (only coupled by the Coulomb force) may influence 
each other. Since the Coulomb repulsion is relatively 
small, a sizeable effect will only arise when the two con- 
ductors are very close to each other. This condition can 
be met in mesoscopic systems — where with the advent 
of new technologies (e.g., carbon nanotubes) the prob- 
lem of Coulomb drag attracted more and more attention 
(e.g., Refs. lll2H 3ll4r) — or in the spin-drag problem. For in- 
stead of considering two conductors one may look at drag 
effects between different fermion species, e.g., fermions 
with different spin orientation. Since fermions with dif- 
ferent spins are not spatially separated, there is a large 
Coulomb force between them which can lead to all kinds 
of correlation effects, e.g., a drag effect. In the last years, 
the interest in spin-dependent transport increased. One 
key problem is the generation of a spin-polarized current, 
i.e., a current where only fermions with one of the two 
spin orientations flow. In this context it is important to 
keep in mind that the spin-polarized current may affect 
the fermions with the opposite spin orientation. Hence, 
the drag effect may play here a crucial role even though it 
is experimentally not directly accessible. (This is because 
the driving potentials are in general not spin dependent.) 



fermions with opposite spin orientation. 

In the Section [H] we present the model and give some 
central definitions for the spin-drag problem. The Sec. 
nil contains the technical details on the subject of the 



For this "spin-drag" problem the trans-resistivity of 
higher dimensional systems has been investigated in pre- 
vious publications, e.g., using the Boltzmann equation 5 
or the random-phase approximation. 6 In this paper we 
focus on the Transconductance for the Hubbard model. 
While most authors used a bosonization approach to 
compute the conductance^ we will use here, for the 
study of the Transconductance, a Monte Carlo method 
which we introduced in a recent paper ^ The strategy we 
followed there was to map our fermion system via the 
Jordan- Wigner transformation to a spin system which 
can be analyzed by efficient though standard Monte Carlo 
techniques. We will now extend this method to the one- 
dimensional Hubbard model concentrating on the ques- 
tion how a spin-polarized current (driven by a voltage 
drop which is assumed to be spin dependent) affects the 
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MC simulations. The MC method of our choie 
a variant of the St ochastic Ser ies Expansion (SSE) as 
introduced in Refs. Illll2ll3ll4 This method allows an 
investigation of the one-dimensional Hubbard model^ 5 *^ 

In the following Sec. IIV Al we present our results for 
the standard Hubbard model and compare with analyt- 
ical predictions from bosonization theory. To obtain a 
spin-polarized current we add in Sec. ITVTSl an impurity 
to the system, which acts like a combination of a one- 
site chemical potential and a one-site magnetic field. We 
show that such a "magnetic" impurity can produce the 
desired spin-polarized current. 

The Hubbard model can also be mapped to a spinless- 
fermion ladder. Hence, our results may also describe that 
situation, but there one might argue that the very spe- 
cific modeling of the Coulomb (on-site) interaction in 
the Hubbard model is unreali stic (and may differ from 
other approaches, e.g., Refs. III EJ3B ). We therefore dis- 
cuss two variants (with additional interaction terms) of 
our model. First, we will discuss in Sec. IIV Cl a situation, 
where fermions with different spins live on different sites. 
The full system has the geometry of a zig-zag chain. Sec- 
ond, we show in the appendix that a spin-polarized inter- 
action leads to equal Cis- and Transconductance. This is 
similar to the "absolute" drag result found, e.g., in Rcf. 
i 



II. DEFINITION OF THE SPIN DRAG IN THE 
HUBBARD MODEL 

Our model Hamiltonian is the standard Hubbard 
model (with N sites or atoms) 

#Hubb = ^ (<4,<t c "+1,<t + cl+l,a c n,o 

-(j,}^n n>a - + U}^n n ^n nt i, (1) 
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where n n ^ is the occupation number of fermions with 

spin a at site n, and Cn, CT is the corresponding annihila- 
tion (creation) operator. To perform our transport cal- 
culations, we will use the approach from Ref. 0- Since 
the hopping term does not connect fermions with differ- 
ent spins, it is natural to consider current and potential 
operators for each spin orientation separately. Following 
Ref. |9| the potential operators read then: (e being the 
charge unit and x being the position of the voltage drop) 

n>x n>x 

The conductance (of a spinless-fermion chain) is the lin- 
ear response of one potential operator to another; there- 
fore the explicit form of the current operators is not 
needed here^ As we have two potentials, we can define 
four transport quantities (conductances) which de- 
scribe the (linear) response of P l to P- 7 where i,j G {T, I 
}. Further details on how to evaluate the g^s are to be 
found in section ITTT1 

For the moment we will discuss only symmetric models, 
i.e., we have spin-rotational invariance (the only asym- 
metric model that we will discuss appears in Sec. IIVB)l ; 
thus, we end up with only two distinct quantities. We call 
9c '■= 9a = fftT the Cis conductance and gt := gii = g-\[ 
the Trans conductance. 

The naming conventions come from the physical inter- 
pretation of these coefficients which is the following: If we 
switch on at a certain time a (supposedly) spin-polarized 
potential which acts only on spin-up fermions (i.e., we 
add a time-dependent perturbation of the form 

vp ] e(t) 

with V being the voltage amplitude and 9 the Heavyside- 
step function) we will find a current of spin-up fermions 

I ] = 9cV. 

This is the drive current governed by the Cisconductance, 
but there will also be a current of spin-down fermions 

I L = 9tV, 

the drag current (governed by the Transconductance). 
The latter may be nonzero, even though the spin-down 
fermions do not feel the applied potential. 

The situation of a nonvanishing Transconductance (or 
drag current) is called Coulomb drag. This problem has 
been studied, e.g., for coupled spinless-fermion systems 
by bosonization in Refs. IJQId an d to second-order per- 
turbation theory in Ref. |j. 

Normally, (since spin-polarized potentials are not 
available) one is only interested in the (full) conductance 
of the Hubbard model, where both fermion species feel 
the same potential, and the full current is the sum of 
the currents of the spin- up and spin-down fermions. As 
is straightforward to see, the Cis- and Transconductance 
give us directly the conductance of the Hubbard model 
via the relation g H ubbard = 2(g c + g t ). 



III. DESCRIPTION OF THE MONTE CARLO 
METHOD 

Before starting with Monte Carlo we have to cast our 
Hamiltonian in a convenient form. 

Using the Jordan- Wigner transform the Hubbard 
model can be mapped to a spin ladder. To each oc- 
cupation operator we introduce a spin operator, i.e., we 
replace n n ^ — > S$ n + 1/2 and n„j — > Sfn+i + 1/2- If wc 
express the Hubbard Hamiltonian with those new opera- 
tors we obtain the following spin ladder (with 2 N sites) : 

n 

+ \P ^2n^2n+l + U Sf n+1 S-j,^] , (2) 

n 

where the sites with even number represent spin-up 
fermions and the sites with odd numbers, spin-down 
fermions. (For the Hubbard model one has to put 
J z = = U'. These parameters are used to model the 
spin-polarized interaction from the appendix and the zig- 
zag chain from Sec. IIV CI see bottom half of Fig. ^ The 
sites 2n and 2n + 1 in the ladder represent therefore one 
atom of the Hubbard model and interact via an Ising 
interaction representing the Coulomb force (see Fig. QJ. 
The hopping amplitudes satisfy J x = 2i, and the strength 
of the magnetic field B in Eq. J5J) is obtained from the 
chemical potential via the relation B = U/2 + U' /2 — /i. 
Half filling corresponds therefore to B = 0. The two po- 
tential operators Pj'^ introduced in the previous section 
can be related to potential operators for the two chains 
(e being the charge unit) 

Px — e J]] £>2n+l> Px — e £>2n- 

n>x n>x 

Then we may obtain the four conductances gij introduced 
in the previous section by computing 6 { J., j}) 

9u(um) = -wk/S I cos{uj M T){P l x P^{iT))dT (3) 
Jo 

at the Matsubara frequencies u>m = 2irM((3H)~ 1 , M g 
N, and then extrapolating to w = 0. (The extrapo- 
lated value should not depend on x or y& We chose 
x = N/2 = y — 1.) For the extrapolation from g{u>M) 
to g(oj = 0) we will use a quadratic fit from the first 
three Matsubara frequencies. [We will use open (OBC's) 
instead of periodic boundary conditions (PBC's).] Since 
the Hamiltonian contains Heisenberg-like interactions as 
well as Ising-like interactions, it is advantageous to use 
the Stochastic Cluster Series Expansion (SCSE) intro- 
duced in Ref. EB For the Hubbard model the SCSE 
gives essentially the same update scheme as the one used 
in Ref. We will explain it now shortly. 

Following Ref. Kj we split the Hamiltonian according to 
H = h, but this time into four-sites clusters, called 
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FIG. 1: The model Hamiltonian that will be discussed in this 
paper. Solid lines indicate a full Heisenberg-like interaction 
between the sites; dashed lines stand for sites coupled only 
by a z-z term (Ising-like interaction). The upper part is for 
U' = 0. The lower for U' = U. (Other values of U' are 
not considered.) A plaquette as used in the MC scheme is 
indicated by boldface lines. 

plaquettes (see Fig. [TJ. The following operators belong 
to the plaquette V (containing interactions between the 
sites 2n, 2n + 1, 2n + 2, and 2n + 3): 

= Jx^2n^2n+2l^ ^"P ~ J x & '2n^2n+2 1 ^ 

i,(3) j c+ a- /o /,( 4 ) _ j q- q+ /o 

= C + U/2S% n S% n+1 +U/2S% n+2 S% n+3 + 

I TT'qz q z i 7 P qz + T q z q z 

' u °2n+l°2n+2 ' J z °2n+l J 2n+3 ' J z°2n°2n+2 

+ B/2(S'2„ + S^n+i + S'2n+2 + ^2n+3J- 



considerably. 

IV. NUMERICAL RESULTS 
A. Transconductance in the Hubbard model 

1. Comparison with the Hubbard model 

In bosonization theory the Hubbard model is described 
by two boson fields <E>f ^ representing the degrees of free- 
dom of different spin orientations. The current operators 
for the spin sectors are then given by Jj^ oc d x §^^m 
The conductance can be written in the form of a current- 
current correlator and may be evaluated in terms of the 
Luttinger-liquid parameters K Pi<T of the charge and spin 
field& <S> p/a = ($ T ± <f>i)/V2. The result is (usin g th e 
linearity of the correlator and the results from Ref. fl7f) 

9c = \{K P + K a ) gt =\{K P - K a ). (4) 

When fi = U/2 we are at half filling, where Umklapp 
processes are responsible for a gap in the system^ The 
(charge) gap A(U) depends on the Hubbard repulsion U 
and is finite for all U. 



2. Numerical simulations 



The set f) consists of all h$ for all plaquettes V and all 
t=l,...,5. 

The heart of the SCSE program is the so called loop 
update, where a spin flip of a subset (loop) of all spin 
variables is proposed. Since the sites with even and odd 
numbers form two chains, which are only coupled by a 
z-z interaction term, we find that the set of spin variables 
that will be flipped in the loop update belongs entirely to 
one of the chains. Therefore, we can view the new algo- 
rithm as making loop updates for each chain separately. 
During a loop update for one chain the spin variables of 
the other chain remain fixed. The consequence of this is 
that, if we update, e.g., the even chain, then operators 
with superscript i = 3, 4 can be neglected (are irrelevant 
for the loop construction), and the coupling terms (be- 
tween the chains) reduce to magnetic field terms (for the 
even chain). 

It is however advantageous to consider another vari- 
ant of the loop update. The construction is similar to 
the first variant, but now we propose spin flips for both 
chains, i.e., the spin variables belonging to sites 2n and 
2n + l are flipped simultaneously. This may be viewed as 
a construction of two parallel loops — one for each chain. 
Since the two loops must be parallel, the number of pos- 
sible transitions between different plaquette states is re- 
duced. This may lead to a less efficient algorithm^ but 
one should note that this parallel-loop update becomes 
deterministic for the case of B = and hence enhances 
the efficiency of the algorithm (at least for this situation) 



We present now Monte Carlo results for the Transcon- 
ductance in the Hubbard model Eq. [or Eq. @ for 
U' = 0]. 

We performed simulations for two different chemical 
potentials: First, [i = U/2 corresponding to half filling 
and second /i = 0. In the latter case the system is no 
longer at half filling, but has a U -dependent filling, which 
is shown in Fig. 1111 of the appendix (in the large-[/ limit 
the system reaches quarter filling). 

We show g c and g t as a function of U for the two dif- 
ferent /i's in Fig. |2 The figure shows that the Coulomb 
drag is very sensitive to a change in chemical potential. 

Let us first look at the half-filled case (Fig. |2J). If U 
is very large the Coulomb repulsion acts as an effective 
projection to those configurations satisfying = —P\ 
This implies that g c + gt — > as U — > oo. This con- 
templation is in accordance with Fig. [21 We should 
actually expect from Eq. (0J that for T = we have 
9c + gt = K p = because of the charge gap A(Z7) > for 
all U > (cf. Ref. El). This should lead to a discontinu- 
ous jump at U = 0, because without the Coulomb force 
evidently g t = and g c is the conductance of uncou- 
pled spinless-fermion chains from Ref. Il7l Here we em- 
phasize that our method is a finite-temperature method, 
which means that the conductances calculated by us in- 
terpolate smoothly between the values for U — and 
U = oo. The crossover is expected to take place at that 
interaction value Ut which satisfies A(Ut) — ksT. It 
is therefore interesting to see how g c t scale with tem- 
perature. However, our method gives only access to the 
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FIG. 2: Cis- and Transconductance (filled/empty symbols) 
of the Hubbard model as a function of U for two different 
fx's (200 sites, T = 0.02t/fc s , OBC's, 2 • 10 s MC sweeps.) 
The (solid) arrow indicates the Ut for which the charge gap 
(present at half filling) satisfies A(Ut) ~ UbT. (The dotted 
arrow shows Ut/2 where A((7t/2) = fesT/2 for comparison.) 



low-T regime^ such that we will compare here results 
for only two different temperatures, fc^T = O.Oli and 
ksT = 0.02<. In the two subsequent simulations for 
the Cis-/Transconductance we did not find any differ- 
ence at all. This implies only a weak temperature depen- 
dence (at low T) for the interaction Ut which governs the 
crossover. Since A(Z7) is known from analytical results^ 
we can calculate the two crossover interactions — defined 
by A([/ . 02 ) = 0.02t and A(Z7 .oi) = O.Olt— finding that 
Uq.q2 ~ 1.25i and £/o.oi ~ 1.12i do not differ much, as 
expected (they are also both indicated by arrows in Fig. 
0. 

Another important consequence of g c + gt — > is that 
the signs of the Cis- and Transconductance are opposite 
or — in terms of the spin-up and spin-down currents — that 
the induced drag current flows in the opposite direction 
of the drive current^ 

Now we will turn to the spin sector. We have K a = 
1 by spin-rotational invariance of the Hubbard modeli& 
implying [see Eq. lfl|] g c — g t = 1 for all U which is very 
well satisfied by Fig. [21 

Putting the two results for g c ± g t together, we obtain 
g c = Q.5e 2 /h = —g t valid at high U. This large-?/ limit 
of g c may also be computed in second-order perturbation 
theory. In this approximation the Hubbard model can 
be mapped to a Heisenberg model. The operator [on the 
Hilbert space of the original Hamiltonian Eq. QJ] Pj, 
which is effectively equal to — P+, is identified with the 
operator (on the Hilbert space of the effective Heisenberg 

model) P x = En> x T n- [ T n = K,T " "nj I 2 is the 

spin operator for fermions; here denoted by T to avoid 
confusion with the spin operators appearing in Sec. IIIII ] 
Applying the results from Ref. the computation of g c 
reduces then to the computation of the spin conductance 
of the Heisenberg model, which equals one half in units 
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FIG. 3: The model Hamiltonian with an impurity at site N/2 
for N = 6 in the spin ladder representation from the upper 
half of Fig. Q (for U' = 0). The site on which the impurity 
potential acts is encircled. 



In the case of zero \i (again Fig. EI) where the system is 
away from half filling there is no charge gap. Hence, K p 
is finite, and so Eq. tells us that g c + gt — K p does not 
decay with U . [Here we note that g c + gt agrees (within 
error bars) with the values for K p available in Ref. ITsl ] 
We still have K a = 1, which leads to g c — gt — 1 for all 
U (again very well satisfied by the figure). 

Finally, we consider the large- U limit. Inserting 
K p (U = oo) = 0.5 and K a = 1 from Ref. [H in Eq. 
(0J yields g c = 0.75 and g t = -0.25 (units e 2 /h). These 
results are in accordance with the figure. (Note that the 
statistical error increases with U such that we cannot 
compute g(U) for sufficiently high U in order to extract 
the large-Z7 limit accurately.) 



B. Magnetic Impurity 

In this subsection we will study the influence of an 
impurity. The modeling of the impurity follows Ref. 0, 
but for the spin-drag problem it is natural to consider a 
spin-dependent impurity, as we will do here. We extend 
therefore our Hamiltonian in the following way 

H = i?Hubb + -Blmp^AT/2,T ' 

i.e., we introduce an (impurity) potential at exactly one 
central site (which acts only on one spin orientation, see 
Fig.EJ). 

Fig. 0] shows Cisconductances and Transconductances 
as a function of the impurity potential Bi mp at half fill- 
ing. (The example is chosen such that the Transconduc- 
tance in the unperturbed system is relatively large.) The 
conductance of the Heisenberg chain with one impurity, 
which is the large-£7 limit of g c , is given for comparison. 

Although the two Cisconductances, g^ and gn, could 
in principle differ (the model is now asymmetric) they 
do not in the case of half filling — at least within error 
bars. Both Cisconductance and Transconductance go — 
more or less linearly — to zero as the impurity strength 
increases. 

We note that within error bars g c = —gt such that the 
full conductance of the system 

9 = 2(3c + 5t) 
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FIG. 4: The Cis- and Transconductance (filled/empty sym- 
bols) as a function of the impurity strength _Bi mp at half 
filling. The impurity acts on spin-up fermions. Note that 
the Cisconductance for the spin-up fermions g-ff (triangle up) 
might differ from the one of the spin-down fermions gu (tri- 
angle down). The two Transconductances are the same. The 
conductance of an xxx (Heisenberg) chain with one impu- 
rity (which should coincide with the large- U limit of the Cis- 
conductance) is given for comparison, ((i = U /2, U = 2t, 
N = 192 sites, T = 0.02i/fc s .) 




FIG. 5: The Cis- and Transconductance (filled/empty sym- 
bols) as a function of the impurity strength Bi mp in the 
Hubbard model away from half filling. The impurity acts 
on spin-up fermions. Note that the Cisconductance for the 
spin- up fermions (triangle up) differs from the one of the 
spin-down fermions gu (triangle down). The two Transcon- 
ductances are the same, {(i = 0, U = At, N = 192 sites, 
T = 0.02t/fc s .) 



remains zero after insertion of the impurity. Further- 
more, investigations with our method at different tem- 
peratures find no sizeable T-dependence. 

In the case of zero chemical potential, (i = 0, we find a 
splitting of the two Cisconductances (see Fig.[3J). This is 
particularly interesting since this finding implies a spin- 
polarized current. If we assume a (spin-independent) 
driving potential of the form 



(V being the voltage amplitude), then the current is 

i = i ] +i [ , ^ = {9u+9nW I 1 = Gm + 9u)V. 

The average spin of a fermion in the current is therefore 
(using 5iT = g u = g t ) 
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different from zero (see Fig. El- 
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FIG. 6: The negative average spin polarization (-S) of the 
induced current as a function of the impurity strength Bi mp 
in the Hubbard model away from half filling. The impurity 
acts on spin-up fermions. The dashed curve is obtained by 
fitting gii, gu and gt in Fig. with an exponential ansatz 
and substituting these fits into the formula for <S. (fx = 0, 
U = At, N = 192 sites, T = 0.02t/k B .) 

An interesting problem is the question whether for 
(i = the Cisconductance in the pure (spin-down) sec- 
tor <?j| survives or not when we increase -Bi mp to infinity. 
Note that a finite gu would imply a total spin polariza- 
tion, i.e., S = —1/2. The limit Bi mp — > oo of g^ can- 
not be taken directly (because of problems with the MC 
simulation), but here we note that there is another way 
to model the impurity. Instead of applying a local mag- 
netic field on one site, one can introduce a weak link, i.e., 
decrease the hopping amplitude for spin-up fermions be- 
tween the sites N/2 and N/2 + 1 from the initial value t to 
tj mp . These two variants of impurities behave similarly^ 



V{P ] + P l ) 



We computed the Cisconductance for the unaffected spin 
orientation for the model with i j mp = (corresponds to 
B Imp = oo) at (i = 0, U = At, N = 200 sites. We find a 
value of about [0.79 ± 0m]e 2 /h. 

The different behavior of the Cisconductance in the 
(unaffected) spin-down sector at half filling and away 
from half filling may be explained as follows: Suppose 
Simp and U are large. The effect of the -Bi mp term on 
the fermions is that it forbids occupation of the impu- 
rity site for one of the two fermion species (in our case 
spin-up fermions). At half filling a spin- up fermion can 
hop only from one site to another by exchanging the 
site with a spin-down fermion (there are no empty sites), 
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i.e., simultaneously with the spin-up a spin-down fermion 
must hop in the opposite direction et vice versa (imply- 
ing g c = —gt). A fermion of a certain spin index can then 
only pass the impurity site if accompanied by a fermion 
of opposite spin (which moves in the opposite direction) . 
Since the impurity site is forbidden for one of the two 
fermion species, no fermion can pass the impurity site, 
and both Cisconductances must go to zero. 

Away from half filling the hopping of an spin-up 
fermion does not necessarily require the hopping of a 
spin-down fermion (the spin-up fermion can hop to an 
empty site) and hence the impurity affects only one of 
the two Cisconductances. 




C. Zig-zag chain 

So far we have dealt with a system of two fermion 
species, where the two species reside on the same set 
of sites. In contrast to this, the bosonization ap- 
proaches considered mostly systems of two coupled 
spinless-fermion conductors. We can compare our results 
with that situation, if we interpret the Hubbard model 
as a spinless-fermion ladder. Here one assumes that each 
fermion species lives on a different conductor [i.e., the two 
indices (n, |) and (n, J.) are supposed to label (spatially) 
different sites; compare Sec. IIIII and upper half of Fig. 

But one should note that for this case the parameter 
U should be small (since the distance between separated 
conductors is large) and the Coulomb interaction should 
be long-ranged (not on-site as in the Hubbard model). 
Hence we are led to the question how a variation in the 
interaction term modifies our results. 

To address this question we add a new interaction term 
to the Hubbard model 

H = i?Hubb + U ^ ^n,4.«n+l,t • 
n 

This Hamiltonian corresponds to Eq. (J2J with U' = U. 
One can justify introducing this new interaction term if 
the fermion species live on different sites where each spin- 
down site (n, I) lies between two spin-up sites, (n, f ) and 
(n + 1, T)- This model has therefore the geometry of a 
frustrated zig-zag chain as depicted in the lower half of 

Fig.n 

One should note that this system has a total of 2N 
sites, N sites for each fermion species. [Although it would 
be useful to adopt the notion of a system with two cou- 
pled (spinless-fermion) chains, we will keep here the no- 
tation of a system of spinfull fermions.] 

The results for the spin drag in this model are shown 
in Fig. [7] We discuss again two chemical potentials: one 
is [i = U implying half filling, the other is again fj, = 0. In 
the latter case the (mean) occupation number per site n p 
is different from one half (the occupation at half filling) 
and depends on U. It is shown in Fig. |H] 

One sees in Fig.[7|that \g c ,t\ grow with the strength of 
the interaction. This may be explained as follows: First, 



FIG. 7: Cis- and Transconductance (filled/empty symbols) of 
the zig-zag chain for two magnetic fields (120 sites per fermion 
species, T = 0.02t/k B , U' = U in Eq. Q, OBC's, 2 ■ 10 5 MC 
sweeps.) 

the Coulomb interaction mediates an attractive nearest- 
neighbor interaction for fermions with equal spin orienta- 
tion (this is a consequence of the frustration). Therefore, 
in a simple approximation the only effect of the Coulomb 
interaction is to renormalize the Luttinger-liquid param- 
eters for the two spin sectors K^'K Since the Luttinger 
parameter for a spinless-fermion chain increases with the 
strength of the attraction^ we expect that IO>* increases 
as U increases. Since K*'^ gives the conductance of one 
spin sector— (which is essentially the Cisconductance) 
we have that g c increases with U. 

One may also infer from the figure that the dependence 
of g c .t on a chemical-potential shift is weak. Within error 
bars g c decreases only slightly upon shifting fi away from 
half filling. 

One should note that in the limit U = oo the ground 
state is a spin-polarized configuration (see Fig. [§} . For 
fi = U this means that all conductances are zero in the 
large-£7 limit, for one spin sector is empty and the other, 
completely filled. In contrast to this, for fi = one of the 
two spin sectors may remain conducting. The crossover 
to the ordered state occurs at values of U larger than 
3t which may be seen by simulating and comparing the 
occupation number for different states. For \i = the 
difference in occupation (between the two spin sectors) 
n <y = (I J2n ( n »,T _ n n,l) |)/(2AT) is shown in Fig. for 
fi = U it is zero within error bars as long as U < 3t. We 
conclude that for the values of U considered in Fig. the 
two spin sectors have approximately the same filling. 



V. CONCLUSION 

In this paper we discussed the spin drag for the Hub- 
bard model at zero temperature. We found that the 
Transconductance is negative — at half filling the Umk- 
lapp even enforces g c = —gt- In that respect our situation 
is different from two coupled Tomonaga-Luttinger mod- 
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FIG. 8: Occupation per site (n p ) and difference between the 
occupations of the two fermion species (n CT ) for the zig-zag 
chain. (100 sites per fermion species, PBC's, fi = 0, 10 4 MC 
sweeps, T = O.lt.) 

(KKK) 

FIG. 9: One of two possible ground state configurations for 
the zig-zag Hamiltonian when U » t. The sites in the lower 
row are occupied by spin-up fermions, the sites in the upper 
row are empty. In the other ground state, the lower row 
is empty, and the upper completely filled with (spin-down) 
fermions. 

els as considered in Refs. which do not incorporate 
Umklapp. The "absolute-drag" result of the form g c = g t 
(e.g., from Ref . 0) can only be recovered by introducing 
a spin-polarized interaction (see appendix). 

If we assume that a given potential is in general not 
spin dependent, the only relevant quantity is the full con- 
ductance g = 2{g c + gt), which is only nonzero away 
from half filling. Here both spin orientations contribute 
equally to the current. However, the situation changes 
when we add a magnetic impurity. Even if the driving 
potential is still spin independent, the resulting current 
will be (partially) spin polarized, if we are away from 
half filling. In the limit of a large impurity potential the 
current will be fully spin polarized. 



APPENDIX A: SPIN-DEPENDENT 
INTERACTION— BROKEN SPIN-ROTATIONAL 
INVARIANCE 

In the Hubbard model Cis- and Transconductance have 
opposite sign, in sharp contrast to the bosonization re- 
sults (for coupled spinless- fermion chains), where Cis- 
and Transcoductance are both positive. The discrepancy 
may come from the different modeling of the interaction. 
In the bosonization approaches each chain is given by an 



interacting system (i.e, the Luttinger-liquid parameter K 
may differ from one) , in the Hubbard model each spin sec- 
tor alone is represented by a noninteracting fermion sys- 
tem. We will show in this appendix that a spin-polarized 
interaction leads to a positive Transconductance as found 
in the bosonization approaches. To this end we will now 
discuss the following variant of our Hamiltonian: 

H = i?Hubb + J ^ n ^ _ l/2)(n„ +1 , CT - 1/2). (Al) 

Here the new J z term breaks the spin-rotational invari- 
ance. Hence K a may now be different from one. 

First we consider the large-U limit at half filling. The 
U term acts then as an effective projection to the configu- 
rations with exactly (because of half filling) one fermion 
per site, i.e., n n ^ = 1 — n n ^. We now set-up an ef- 
fective (second-order perturbation theory) Hamiltonian. 
From the kinetic-energy term we get again a Heisenberg 
model with exchange parameter At 2 /U. The J z term of 
the Hamiltonian does not change the configuration (in 
the occupation-number basis) and gives therefore a di- 
rect energy contribution 2 J z J2 n T*T* +1 to the effective 
Hamiltonian where T£ — (n n ^ — n nt {)/2 denotes the spin 
of the fermion on site n. The full effective Hamiltonian 
reads 

H cS = ]T(4i 2 •■/U)(T+T- +1 +T+ +1 T-)/2 

n 

+ ^(2J z +4t 2 /c/)r^ +1 , 

n 

and is an xxz chain. If the anisotropy is larger than the 
hopping amplitude, i.e., if J z > 0, this model is gapped 
(implying both a charge and a spin gap in the original 
model) . We therefore expect that Cis- and Transconduc- 
tance go to zero, if we increase U and keep a finite J z . 

Now we consider a zero chemical potential /i — 0. We 
expect that this chemical potential shift away from half 
filling closes the charge gap, but leaves the spin gap more 
or less unaffected. We consider again the large-U limit. 
In any configuration the J z term of the Hamiltonian gives 
the following contribution for two neighboring sites 

— J z /2 if the two sites are occupied with anti-parallel spins, 

J z /2 if the two sites are occupied with parallel spins or 
are both empty, 

if one site is occupied and the other, not. 

We assume that there is a spin gap and that the (degen- 
erate) ground state configurations are those for which the 
spins of the particles are ordered antiferromagnetically. 
If only these configurations are allowed, the J z term can 
be represented as a one-site potential with a contribution 
±J 2 /2 for empty /occupied sites. (One obtains the same 
energy contributions as from the J z term of the original 
Hamiltonian, if one keeps in mind that each site appears 
in precisely two pairs of neighboring sites.) We can set- 
up the following effective Hamiltonian (this is just the 
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FIG. 10: The same as Fig. H but for J z = 0.8t (here N = 
140). The simulations away from half filling suffer from large 
autocorrelation times. 



Transconductance in Fig. [23 In the large-?/ limit we 
find good agreement with our prediction that g c — 9t — 
0.25 e 2 /h which gives credit to the simulation data de- 
spite the large autocorrelation times. 

Here we want to stress once again the remarkable fact 
that the sign of the Transconductance (the direction of 
the induced current) changes when we switch the mag- 
netic field and the spin-polarized interaction on. (The 
Transconductance is for all U negative in Fig. [21 whereas 
in the present situation we expect gt = g c = K p /2 > 
for T = 0, U = oo.) 

Occupation in the ground state — Since we identified the 
ground state of the Hamiltonian Eq. (|A1(I H([i = 0,U — > 
oo) with the ground state of the xx chain in magnetic 
field, we can calculate the occupation per state of this 
Hamiltonian in the large [/-limit, the result being: 



restriction of the original Hamiltonian to the assumed 
ground state configurations, i.e., zeroth order in U) 



E 



[(R+R 



n+l 



\n n>a )/N = 1 - arccos(J z /[2i])/7r. 



This prediction may be tested against a Monte Carlo 
simulation. We find good agreement (see Fig. ITT|) . 



which is an xx chain in magnetic field, where the "spin" 
operator R* (this time denoted by R to avoid confusion 
with previous spin operators) is defined by R z n = —1/2 
if there is a particle on site n and R^ = +1/2 if site 
n is empty. Since the effective Hamiltonian describes 
the charge part of the Hamiltonian the full conductance 
2(g c + gt) for the original model should coincide with the 
conductance of the new Hamiltonian which is e 2 /h as 
the system is noninteracting£i£ Since we have g c = g t 
by the assumption of a spin gap and Eq. , the relation 
g = 2(g c + g t ) yields g c = 0.25 e 2 /h = g t . 

In principle the model Eq. (|A1|) can also be analyzed 
with the Monte Carlo method developed in this paper, 
but we found that the simulation for this case is problem- 
atic: We measured large autocorrelation times for finite 
J z and \x (e.g., for the computation of the compressibil- 
ity). We therefore must restrict ourselves to J z < 0.8t. 

For J z = 0.8t we present results for the Cis- and 




FIG. 11: Occupation per state (away from half filling) for 
the Hamiltonian Eq. 1A1I H(fi = 0) for different J z . The 
predicted high U values are given as dashed lines. (500 sites, 
PBC's, 10 5 MC sweeps, T = O.li.) 
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